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Abstract 

In this paper we construct multiparametric families of two dimensional metrics 
with polynomial first integral. Such integrable geodesic flows are described by solu¬ 
tions of some semi-Hamiltonian hydrodynamic type system. We find infinitely many 
conservation laws and commuting flows for this system. This procedure allows us 
to present infinitely many particular metrics by the generalized hodograph method. 
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1 Introduction 

The problem of integration of geodesic flows on a two-dimensional surface appeared in 
classical mechanics very early and was extensively investigated in the XlX-th century 
focusing mostly on local approach. The XX-th century was more centered on global 
behavior, problems of local description of trajectories were not that much popular. An 
overview of both aspects may be found in numerous sources, we give here only two ref¬ 
erences [2, 3]. Partially this loss of interest to the local problem was not only due to the 
importance of the global problems; it seems reasonable to ascribe this loss of interest to 
absence of new ideas for the local integrability problem. The situation, in our opinion, 
has changed in the last decades after publications [5], [1] where the authors had remarked 
that the equations for the coefficients of first integrals polynomial in momenta for two 
specific low-dimensional cases belong to the class of diagonalizable hydrodynamic type 
systems integrable by differential-geometric means; the appropriate theory for such non¬ 
linear systems of PDFs was developed in the very end of the XX-th century (cf. [4, 14]). 
In [12], developing the preliminary results of [5], we demonstrated how to apply the tech¬ 
niques of integrable hydrodynamic type systems to the case of so called one-and-a-half- 
dimensional systems (one-dimensional mechanical systems with the potential depending 
on time). Below we investigate (using a bit more sophisticated technologies) the problem 
of local description of two-dimensional Riemannian metrics with geodesic flows possessing 
a polynomial first integral of arbitrary high degree. 
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In [1] the N component hydrodynamic type system 

a° = af = + [{k + — (A^ + 1 — (1) 

where k = 1,.. ., iV — 1 and = 1, was derived as the system for the coefficients a^{x, t) 
of a (global smooth) polynomial hrst integral^ 

( 1 ^ ^ 

F[x, t, pi, Pa) = ^ N-k Pi~^P2, (2) 

k=0 ^ 

for metric in semi-geodesic coordinates 

ds'^ = g‘^{x, t)dt^ + dx"^ (3) 

on a 2-dimensional torus, where g = a^~^. As shown in [1], any smooth metric on a 
2-dimensional torus possessing a polynomial first integral can be reduced to such a form. 
The metric (3) corresponds to the Hamiltonian H{x,t,pi,p 2 ) = \{p\/g‘^{x,t) +P 2 ); the 
system (1) is equivalent to 

dF dH dF dH dF dH dF dH _ 

’ dpi dt dt dpi dp2 dx dx dp2 

The theorem proved in [1] states that the hydrodynamic type system (1) is diagonalizable 
(i.e. possesses a complete set of Riemann invariants) and semi-Hamiltonian so may be 
integrated by the Generalized Hodograph Method (cf. [14]). However the authors of [1] 
did not give a constructive description of the necessary complete set of hydrodynamic 
conservation laws and commuting flows for (1). 

In this paper we construct N inhnite series of conservation laws and commuting flows. 
Thus one can construct a rich multiparametric family of particular solutions to (1) by 
the Generalized Hodograph Method as described below in Sections 2, 7. Our interest is 
focused on local properties of the system (1) and the respective coefficient g{x,t) in (3). 

The structure of the paper is as follows. In Section 2 we discuss the semi-Hamiltonian 
property of hydrodynamic type system (1). We construct a generating function of con¬ 
servation laws and hnd the equation of the associated Riemann surface. We rewrite 
hydrodynamic type system (1) in a diagonal form. In Section 3 we rewrite hydrodynamic 
type system (1) via characteristic velocities and derive analogues of the Lowner equations 
and the Gibbons-Tsarev system. Also we remark that in the two component case hydro- 
dynamic type system (1) is nothing but the simplest two-component linearly degenerate 
hydrodynamic type system (which can be written in appropriate field variables in the 
form Ut = vux,vt = uVx), whose general solution can be presented in implicit form with 
explicit dependence on two arbitrary functions of a single variable. Such a system has 
global solutions. This result conforms to the well-known fact that geodesic flows with 
quadratic hrst integrals can be integrated by separation of variables. In Section 4 we 
introduce new held variables and show how the equation of a Riemann surface rewrit¬ 
ten in these held variables 6^ helps in constructing N inhnite series of conservation law 

^The factor (—1)^ was omitted in [1]. 
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densities. In Section 5 we introduce two hydrodynamic chains as inhnite sets of equations 
compatible with hydrodynamic type system (1). This procedure allows us to construct 
inhnitely many conservation laws (the so called Kruskal series) in a compact form. In 
Section 6 we present the way to construct inhnitely many higher commuting hows and 
inhnitely many associated conservation laws. In Section 7 we adopted the Generalized 
Hodograph Method for construction of a rich inhnite-parametric family of particular so¬ 
lutions for the case of hydrodynamic type system (1). Finally in Conclusion (Section 8) 
we discuss the problem of completeness of the constructed series of conservation laws and 
briehy expose further perspectives of integrability of two-dimensional geodesic hows. 


2 Semi-Hamiltonian Systems and their Integration 

Integrability of a diagonalizable hydrodynamic type system (similar to (1)) means that 
such a system possesses N Riemann invariants r = {r^,, r^), has inhnitely many hydro- 
dynamic conservation laws and locally any solution can be constructed by the Generalized 
Hodograph Method: 

Wi{r\...,r^) = Vi{r^,...,r^)-t + x, i = l,...,N, (4) 

where ^^(r) are the characteristic velocities of the system (1) in the diagonal form 

rl = Vi{r)rl, i = l,...,N, (5) 

and Wi{r) are the coefficients of commuting with (1) hows 

rl = Wi{r)rl., i = l,...,N, (6) 

(no summation over repeated indices is assumed anywhere in this paper). 

Such systems were called “semi-Hamiltonian” in [14] where the detailed exposition of 
the respective diherential-geometric theory was given. 

In this Section we present the equation of a Riemann surface A(g, a) associated with 
the hydrodynamic type system (1). Its branch points r* = A|q=q., where qi are solutions of 
the algebraic equation Xq = 0, are the Riemann invariants of (1). Also we briehy remind 
how to hnd a rich inhnite family of conservation laws and commuting hows (6) for (5). 

Theorem [1]: Hydrodynamic type system (1) is semi-Hamiltonian. 

In this Section we give an alternative proof of this Theorem obtaining the associated 
Riemann surface playing the key role in our construction of the explicit formulas for the 
conservation laws and commuting hows for (1). 

Proof: According to [1] hydrodynamic type system (1) can be equally derived for the 
coefficients a^{x,t) of a polynomial hrst integral 

N 

F{x,t,p) = ^(-l)^a^(a:,t)/(l -/)^ (7) 

k=0 

^The corresponding formula in [1] contains misprints. 
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( 8 ) 


of Hamilton’s equations^ 

, dH , dH 

op ox 

where the effective Hamiltonian function is 

H{x,t,p) = —a^~^{x, 


This means that in fact we have a one-dimensional mechanical system with Hamiltonian 
depending explicitly on the “time” t. Such systems were called “1.5-dimensional” in [7]. 
An integrable subclass of 1.5-dimensional systems was studied in [12] where the general 
technology of hydrodynamic reductions of integrable nonlinear hydrodynamic chains was 
used for explicit description of various integrable 1.5-dimensional cases. In the present 
paper we develop a more general approach for the case studied. 

In this 1.5-dimensionaI setting the corresponding Liouville equation^ 


f=u 

* ’ dp dx dx dp 

takes the form 

= + + (9) 

where instead of the variable p we introduce the auxiliary variable q connected with p via 
the point transformation 


q = 


P 


p = 


q 

V^ + q^ 


( 10 ) 


Substitution of a general ansatz f{x,t,p) = \{q,a{x,t)) with a set of some (formally un- 
specihed) held variables a = {a^{x,t),..., a^~^{x, t)) into (9) leads to an overdetermined 
compatible systems of equations on dX/da\ dX/dq if a^{x,t) are solutions of hydrody¬ 
namic type system (1). Now X{q,a.{x,t)) can be found explicitly: 


A(g,a) 


N-l 


(1 + g" + 


*a‘ 


A:=0 


( 11 ) 


In fact (11) can be obtained by a direct substitution of (10) into (7). Thus a Riemann 
surface with parameters (A, q) used below for the explicit constructions of the conserva¬ 
tions laws and commuting hows is dehned. And vice versa, substitution of (11) into (9) 
and splitting w.r.t. q leads to (1). In a generic case the algebraic equation Ag = 0 for A 
of the form (11), i.e. 

N-l 

[(IV - = Nq^-\ (12) 

k=0 

^This means that geodesic flows can be written either as stationary Hamilton’s equations with two 
degree of freedom or alternatively as non-stationary Hamilton’s equations with one-and-a-half degree of 
freedom. Here the “prime” means the derivative with respect to “t”. 

^The Poisson bracket is standard {/, H} = fpR^ — fxHp- We remind that f{x,t,p) is a first integral 
for the given Hamilton’s equations, for example f = F. 
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has N simple distinct roots qi = gi(a). Since Ag = 0 in these points qi = qi{a{x,t)), 
Liouville equation (9) leads to the hydrodynamic type system (1), written in the diagonal 
form 

rl = (13) 

where the Riemann invariants r* = X\q=qi are determined by (11), i.e. 


N-l 


R(a) = X{qi, a) = (1 + (qi^) (g^)^ + ^(g*) 


mm 

(Jj 


m=0 


(14) 


It is easy to see that the constructed r*(a) are functionally independent since the character¬ 
istic velocities in (13) are distinct. (Functional independence of the velocities Vi = a^~^qi 
will be shown below in Section 3.) 

Under the functional inversion transformation A = A(g, a) —)■ g = g(A, a) the linear 
equation (9) transforms to 

qt = a^~^qqx - (1 + (15) 

which is equivalent to the equation^ 

Pt = (^a^"Vl (16) 

for the so-called generating function of conservation law densities p = p{X, a°, a^, ..., a^~^) 
where A is a parameter. Inhnitely many conservation laws can be found directly from the 
equation of Riemann surface (11) expanding the inverse function g(A, a) (as g —?■ oo, while 
A —1) w.r.t. the powers of (A — 1) and substituting into the second relationship of (10) 
and then into (16).® Since coefficients of the resulting expansion of p{X, a) are conservation 
law densities (see (16)), and the hydrodynamic type system (1) is diagonalizable (see (13)), 
we conclude that (1) is semi-Hamiltonian. The Theorem is proved. 

The next important step in integration of semi-Hamiltonian system (1) is construction 
of the necessary complete set of Wi{r) (depending on N functions of one variable, cf. 
[14]) in order to be able to construct any solution of (1) locally using the Generalized 
Hodograph formula (4). This problem will be discussed in Sections 6, 7. In the following 
three Sections we consider different forms of the system (1) and the associated Riemann 
surface (11) in more detail in order to expose the necessary techniques. 


3 More on Characteristic Velocities and Riemann In¬ 
variants. Triviality of the Quadratic Case (N = 2 ). 

Since the held variables in (1) are connected with the characteristic velocities (the 
eigenvalues of the system matrix) algebraically (see (12) and (13)), one could try to 

®Under the potential substitution p = Sx this equation reduces to the Hamilton-Jacobi equation 
St = — (Sx)'^ for the mechanical system with the specified Hamiltonian H{x,t,p). 

®This remarkable class of conservation laws is usually called Kruskal series. We describe its construc¬ 
tion in more detail below in Section 5. 
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find explicitly the dependence of the functions on the characteristic velocities ni(a) = 
a^“^gj(a), if nj(a) are functionally independent. Thus automatically the hydrodynamic 
type system studied in this paper would be written in a symmetric form. 

Indeed, (15) yields directly the same hydrodynamic type system (1) written as 

= + (17) 


where we just replaced q by the corresponding roots 

Remark; This replacement (g —)■ g^) is not so obvious, because roots g^ of algebraic 
equation determined by the condition Ag = 0 do not depend on the parameter A, while 
the function g depends on it. In fact gfc(a) = g(A, a)|;,=rfc(a). For this reason we give here a 
precise derivation of the above equations. The general theory of hydrodynamic reductions 
developed in [6] implies that g(A,a(r)) satisfies the overdetermined system (the so called 
Towner type equation) 

d^q = d, = d/dr\ (18) 

Qi-q 

where we denoted A = — lna^“^. These equations are in fact the equation (15) for 
g(A,a(r)) after substitution of (13) and splitting w.r.t. rj,. Compatibility conditions 
dk{diq) = di{dkq) for (18) yield the so called Gibbons-Tsarev type system 

d,qk = ^^d,A, d,dkA = 2 }"^'^% d,AdkA, i^k. (19) 

qi - qk [qi - qkY 


One can immediately verify that if g^ are determined by (19), then (17) is fulhlled. This 
computation shows that the hrst (left) part of the Gibbons-Tsarev type system can be 
formally obtained from (18) by replacement g —?■ g^ as well as (17) from (15). 

Taking into account the relationship (see (13)) between characteristic velocities Vk = 
o>^~^qk and roots g^, the hydrodynamic type system (17) takes the form 


{vk)t = Vkvl + Vka, 


N-2 2(nfc)^ + (2a^ ^ - N)vk + (a^ 

X nN — l ’ 


( 20 ) 


where ^ and a" ^ can be found from (12), written in the form 


N-l 


k=0 


5^ (iv - fc) 


(V, 


\fc+l 


\fc-l 


(QAf-l)fc+l 


^N-l\k-l 


\N-l 


= N 


N-1\N-1 ■ 




j = This is a linear algebraic system on unknown functions a^(v), k = 

0,..., {N — 2), except the latest coefficient a^”^(v). One can find, for instance, that 


a 


N-2 


N 

~2 



771=1 


while ^(v) satishes different algebraic equations for different N. For instance, if iV = 2, 
then {a^Y = —Vi, if = 3, then 
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if = 4, then 


(aY+3^2(a3)' + ^4 = 0; 


if = 5, then 


,4^4 




8K 


3(5-i/i)' ' 3(5-i/i) 


= 0 , 


etc. Here Vk are the elementary symmetric polynomials of vj of degree k found by the 
Vieta Theorem from the formal equation 


yk - + V2v’^-^ - Vsv'^-^ + ... = 0, 


so Vi = J2^m,V2 = |[(E^m)^ - ••• 

Taking into account (14) we see that all Riemann invariants can be explicitly 
expressed via the characteristic velocities Vm- If N = 2, we get 


r 


1 




( 21 ) 


where the corresponding equation of the Riemann surface becomes 


A(g,v) 






V^-V2 


if A^ = 3, then 


= 


(Vkf + S^iVkY + VsVk + (V2 + 




3/2 


if A^ = 4, then 


R = 1 + 


(a=^)'(^A:)=' - IViiaYYvkY - V,Vk + - \{aY^V, + iRR 

{{vkY + (a3)2)2 


However, inverse formulas for ^^(r) are much more complicated. Only if A^ = 2, the 
inversion is simple (see (21)): 

V2 = 2 - 2r\ ui = 2 - 2r^. 

Thus, hydrodynamic type system (1) in the two component case 

a° = a^al, a] = aR° + 2(1 - a°)a^ 

is linearly degenerate: 


where 


ri = 2{1 - r^y„ r? = 2(1-r‘)rj, 

— 1, (a/^ — —4(r^ — l)(r^ — 1). 


( 22 ) 

( 23 ) 











Remark: Algebraic equation (12) determines N roots qk{si). However, one can see 
that inverse expressions a^(q) cannot be found if N is even due to degeneracy of the map¬ 
ping (a°,..., —> (gi,..., q^) for even N. But if N is odd, the inverse expressions 

a^(q) can be found. For instance, if = 3, 


a 


1 


3Q3 0 _ Q2 + 2 2 _ 3 

Qi + 2Q3 ’ Qi + 2Q3 ’ Qi + 2Q3 


where Qi = qi + q 2 + q 3 ,Q 2 = gig 2 + gig 3 + <? 2 g 3 , Qs = qiq 2 q 3 - The corresponding expression 
for the equation of the Riemann surface (11) is 


A(g,q) 


(1 q^) 

Qi -\- 2 Q 3 


[(Qi + ‘2Q'i)q^ + 3g^ -l- ^Qd,q + Q 2 + 2]. 


3.1 Linearly Degenerate Case {N = 2) 

Two component hydrodynamic type system (22) corresponding to quadratic hrst integrals 
(2) of geodesics on a two-dimensional surface is linearly degenerate (see (23)), that is 
dvi/dr^ = 0 and dv 2 /dr‘^ = 0. Its general solution r^{x,t), r'^{x,t) depending on two 
arbitrary functions f3{u) and 7 (u) of one variable is well known (cf., for example [13]) and 
may be presented in implicit form: 

t = (3'{u) + 7 '(u), X = (3{u) - u(3'{u) + ■j{v) - vy\v), 

where (for simplicity) we denoted u = 2(1 — r^),v = 2(1 — r^). 

The fact of complete integrability in the case of quadratic hrst integrals is also well 
known, its exposition can be found for example in [2] and [3, Ch. 11], Actually the 
standard exposition of this case (for the isothermic form ds"^ = (/(u) -l- g{v)){dv? + dv^) 
of the metric) in the references given above is equivalent to our result for the metric in 
semi-geodesic coordinates (3). 


4 Symmetric Field Variables and N Principal Series 
of Conservation Laws 

In the previous Section the hydrodynamic type system (1) was written in symmetric form 
(20). However such a symmetric form is non unique. Let us introduce the roots of the 
polynomial 

N-l N 

7 + = JJ(5 _ 6*), 

k=0 k=l 

SO all held variables become elementary symmetric polynomials of new held variables 
6"'. For instance, 

N 

aN-1 ^ (24) 

k=l 
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Thus, the equation of the Riemann surface (11) takes the form 


N 

A = (1 + (25) 

k=l 

After substitution of this expression into (9) the hydrodynamic type system (1) reduces 
to another simple symmetric form 


N / ^ \ 

m=l \n=l / 


( 26 ) 


Indeed, hydrodynamic type system (26) can be derived in three steps: 

1. compute the partial derivatives of In A with respect to the independent variables 


N 


N 


Um 


Nq 


N 


m=l 


m=l 


1 + q‘^ 




m=l 


2. substitution of these derivatives into (9) leads to 


^ Um ^ um 

—— = 

q — q — ]j 

=1^ m=l^ 


E 


■‘'E:Eyr+R‘'-E 


q^ + I 

-^^q - 


,7V-1 


or 


^ um N 

= a^- 


N 


N 


N 


q — b' 

m=l 

which simplihes to 


+ a«->E-^ - ( E*” + E-^ ^ 

m=l m=l^ 


\.m=l m 


^ q — b'^‘ 


N-1. 


E 

m=l 


bT + [(b"'? + - a^-^b^b'^^ 


N 


q — h^ 


= a 


N- 


‘E5"-E5^ 

m=l m=l 


m^N—1 

(^X ■ 


3. Taking into account (24) we conclude that the r.h.s. of this equation vanishes so 
splitting the partial fraction decomposition in the l.h.s. of it w.r.t. q we get precisely 
(26). 

The variables b^ are very convenient for explicit computation of commuting flows and 
conservation laws for the system (1). Since we will need in fact not only the conservation 
law densities (for example p appearing in the left hand sides of equations like (16)) but 
also the fluxes — expressions inside the (.. .)x in the right hand sides of such equations 
— we will call conservation laws the equalities (16) themselves. 

Now we are ready to explain our method of deriving explicit formulas for conserva¬ 
tion laws of (1) using the Riemann surface associated to this integrable hydrodynamic 
type system. One should note that N component semi-Hamiltonian hydrodynamic type 
system has inhnitely many conservation laws and commuting flows — both families are 
parameterized by N arbitrary functions of a single variable (see details in [14]). In many 
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interesting cases this functional dependence cannot be presented in explicit form. Never¬ 
theless N inhnite series of conservation laws and commuting flows can be constructed, for 
instance, if the corresponding equation of associated Riemann surface is known. As one 
can prove (cf. for example [11]) for a vast class of semi-Hamiltonian systems such series 
form the complete basis for the (inhnite-dimensional) linear space of all conservation laws 
(or commuting flows). Another techniques for proving completeness can be found in [12]. 

N infinite series of conservation laws for (26) can be found in three steps: 

1. Let us expand^ q with respect to the local parameter A at the vicinity of each root 
b^\ 

g(^)(A) = 6^ + Agf^ + A2gf^ + A3gf^ + ..., A ^ 0. (27) 

All coefficients can be found recursively, for instance 

ik) ^ (1 + 

n ' 

m^k 


2. Substitute the series into (10), obtaining expansion 

(fc) 


where 


(k) ^- ^ ^ =h^ + XpP + AVf ^ + ^¥ 3 ^ + ..., (28) 




= - , (29) 

y/TTW 

ik) 

Then again all conservation law densities pm can be found recursively, for instance 


(1 + ( 6 ‘) 2 )(»- 3)/2 

n - '>”) ' 

m^k 


(30) 


3. Substitute p^^^(A) into the generating function p of conservation laws (16) 

N 


(p«(A)), + L/TAjAwE'’" =° 


n=l 


and expand both sides w.r.t. the parameter A —?■ 0. Matching the coefficients of the 
same powers of A one obtains N inhnite series of conservation laws, namely 


N 


Uk {k) ^ 

(h% + (yi - ) = 0, iXX). = I 


n=l 




(31) 


^This expansion is given by the so called Lagrange-Biirmann inversion (see [12] for an example of its 
use). 


11 













Remark: As a by-product we obtain the hydrodynamic type system (1), (26) in a 
symmetric conservative form 

/ ^ vm \ 

{h\ = j . i = 1..... iv, (32) 

where we utilized the inverse point transformation (cf. (10), (29)) 


yi-(v)"' 

Its hrst N conservation laws are (31) expressed in variables using (29). 

Alongside with the constructed N series of conservation laws one can obtain another 
(incomplete) set of conservation laws usually called Kruskal series. Its form is much 
simpler and is symmetric w.r.t. the variables U. We give Kruskal series below in Section 5. 


5 Hydrodynamic Chains and Kruskal Series of Con¬ 
servation Laws 


If we introduce new variables (called moments) 


Bk 


N 


k + 1 




k+1 


k = 0,1, 


m=l 


(33) 


then hydrodynamic type system (26) implies inhnitely many equations 

= (V + 2R^)R° - Rf = (A;R^-^ + (A: + 2)R^+^)R°-R°R^+\ A: = 1,2,...; 

(34) 

just hrst N of them are independent since due to (33) all higher moments 

are polynomial expressions of the hrst N moments so hydrodynamic type system (26) in 

moments B^ takes the form 


B^ = {N+2B^)Bl-B^Bl, B^ = (A:R^-^ + (fc+2)R^+^)R°-R°R^+\ k = 1,2 ,..., N-2, 

(35) 

Rf-i = [{N - 1)B^-^ + {N + 1)R^(B)]R° - B\B^{B))^, 

where the dependence i?^(B) = B^{B^,... ,B^~^) is be found from (33) using the 
standard combinatorial results on symmetric polynomials. For instance, 

1. if V = 2, then 

R2(B) = - -{By-, 

6 

2. if V = 3, then 

R3(B) = ^ ^ L{By-, 
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3. a N = 4, then 

B\B) = + B^B"^ - ]^{B^fB‘^ - ]^B^{B^f + ^{B^B^ - ••• 

Nevertheless, one can consider infinitely many equations (34) without above restric¬ 
tions. Similar infinite chains of quasilinear first-order equations (called hydrodynamic 
chains) are very useful for integration of various semi-Hamiltonian systems appearing in 
applications. An overview of this approach can be found in [10, 12]. 

The associated Riemann surface (25) can be expanded at infinity as g —)■ cxo and A ^ 1. 
For convenience we replace below A by p = — In A, so /i —)■ 0 and 

N ^ 

= y ln(l + g^) - ln(g - b^). 

k=l 


Then (see (33)) asymptotically 


N 


= + A +E 


Bk 

nk+l 


k=0 


E 

k=0 


Ck 
nk+1 ' 


g —OO, 


(36) 


and introducing new field variables using their relation to B^ in (36) we obtain another 
remarkable hydrodynamic chain 


Cf = {kC^-^ + {k + 2)C''+^)C° - C^C^+\ fc = 0,1,2,..., 


(37) 


where = B^^ and = B 


2k-l 


\k N 


Embedding of hydrodynamic type system (26) into the hydrodynamic chain (37) for 
arbitrary N allows us to find Kruskal conservation laws® in a compact form. First we 
(using the Lagrange-Biirmann inversion) get from (36) the asymptotic decomposition of 
g and p as /i —)■ 0: 


(JQ Ql 


(C‘)" 


{O 


0'|2 


p(p) = -1 + 




(CO)® 


+ /i" 






0'i3 


(38) 


(C 


OU 


(C' 


0'i5 


2{C 


0\2 






OM 


+ /i" 


^ 5(C1)2 


+ .... 


(39) 


((^0)5 2(C'0)6 8{C^)\ 

Substitute the last expansion into (16), note that = —C^ and equate the coefficients 
at equal powers of jj, obtaining the Kruskal series of conservation laws. The first few of 
them are: 


((C»)-^), = l2C‘(C»)-^]„ 1C‘(C»)-‘|, = ([^ 

5(C1)2 3 


OM 






(^0)5 2(C0)6 ' 8(C0)4 

®We call this asymptotic expansion Kruskal, because M. Kruskal was the first who introduced a similar 
construction for the KdV equation. 


3C1 


(CO)® 2(C0)2 

bic^f 


(40) 


2(C" 


OM 




0'|5 


{C 


0'|6 




0'|4 
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6 Commuting Flows 

The most important part of the integration procednre for a semi-Hamiltonian system (as 
described briefly in Section 2) is construction (preferably in an explicit form) of sufficiently 
many commuting flows, either in a diagonal form (6) or in non-diagonal representation in 
terms of the variables or 6^. We choose the latter possibility in this Section; the next 
Section is devoted to a suitable modification of the Generalized Hodograph formula (4) 
to the non-diagonal form of the commuting flows. 

In the case considered in this paper the conservation laws for the generating function 
have the form (see [8]) 

Pt = (^(p,&^&^...,&^))x 

where the flux S{p, 6^, 6^,, b^) can be represented in a much simpler form S{p, B^) (cf. 
[9]), where i?°(b) is the “zeroth” moment (of the corresponding hydrodynamic chain (34)). 
In this paper we consider the case S = Si{p,B^) = —— p"^. The situation with 
higher commuting flows for the chain (34) is as follows: the first commuting flow has the 
flux (corresponding to the same generating function p{h, A)) of the form S 2 {p, B^, B^), the 
second commuting flow has the flux S' 3 (p, B^, B^, B"^) etc. Assume that all these functions 
Si{p, B^), S 2 {p, B^, B^), S^i^p, B^, B^, B^), ... are the coefficients of an expansion of some 
more complicated function w.r.t. an extra parameter (, we come to the ansatz that such 
a function should be written also in a compact form G{p{\),p{()) (see a more general 
exposition in [10]). 

The conservation laws for the generating function of all higher commuting flows were 
found in [10]. For this particular case they can be represented formally as 


drir,)P{p) = d^G{p{fi),p{r])), 


(41) 


where /i = — In A, p{ri) is obtained by formally replacing the parameter p by the parameter 
7] and 


G{p{p),p{r])) 


V^-p^jp) , 1 . 1 pjy) +1 , . pjp) -pjy) 

^1 -p‘^{ri) 2 p{vi) - 1 -p‘^{p) + -p‘^{p) 


(42) 


The so called “vertex” operator dr{n) is not yet determined and should be specified sepa¬ 
rately for different cases below. 


6.1 The Kruskal Series of Commuting Flows 

For this series infinitely many respective fluxes for the generating function of conservation 
laws of the higher commuting flows can be found substituting the following asymptotic 
expansions for r; —0 (cf. (38) and (39)) 


q{y) 


GO 


icr\ yyl 

{GY) \{GY 


3GY^ 2{GY\ 


+ ..., (43) 


p{y) 


y^ 3 

1 H--h v 

2{GY {CY 


/ G^ 5{GY 
V 5(^ 


s { gy ) ^ 


(44) 
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into (42) and (41) and specifying the expansion of the vertex operator 


= \nr]dto + -Sp + (9p + T]dt3 + + ... 

T] 


to match with the expansion 


G{p{f^),p{v)) =p{p) In?? 


CVl-P^(h) 


P 


+p 




^^/l-p‘^{p) -p{p) lnC'° + -p(/i)ln2 

{cr 


\ _ (ji 

+ ^J - 


c'Vi-p'(h) 


((^ 0)2 ((^ 0)3 

So now we can identify = x and = t, so 

{p{p))to = {p{p))^, (p(/i))p = - (^CVl 

while higher conservation laws are (cf. (16) and (24)) 

(p(p))t 2 = (h) - P{P) In <g° + ^ In ^ I - p(p) In 2^ , 


+ .... 


(p(p))t3 = 


(C'l)^ 


(C0)2 (C'0)3 


\ _ (Jl 

+ ^j (^pm- 


cvi-p"(p) 


Corresponding higher Liouville equations are associated with higher commuting hydro- 
dynamic chains. For instance, the hrst higher Liouville equation 


C 


U = [(1 + + l- \n2] f^ + {q + l)fq [ q\nC^ + — ] 


(JO 


C^\ 


C^J. 


is associated with the hrst higher hydrodynamic chain commuting with (37): 

C% = + ^67^+' - (In67° - 1 + In2) 67^ - [{k + 2)C^+^ + kC^-^]^ 


+ [C\kC^-^ + {k + 2)67^+^) - 67°((A: + 1)67^' + {k + 3)C^+^)] 


6.2 N Principal Series of Commuting Flows 


(JO 


67° 


(C 


0'|2 ■ 


As we remarked above any N component semi-Hamiltonian hydrodynamic type system 
has inhnitely many conservation laws and commuting hows parameterized by N arbitrary 
functions of a single variable (see detail in [14]), but in many interesting cases only N 
inhnite series of conservation laws and commuting hows can be constructed (cf. [11, 12, 
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14]); they usually form a complete basis in the respective linear spaces. Completeness of 
N series constructed here is discussed in Conclusion. 

In this case we again (as in the previous sub-Section) start from (41) and (42). 

The fluxes for the generating function of conservation laws of corresponding higher 
commuting flows can be found from the equations written in a conservative form (cf. 
(32)) 

= 9.G(fti,p(C)). (45) 

Substitution (see (28)) of the expansion of p{() as C —t 0 leads to construction of higher 
commuting flows. However, this derivation is not so trivial. Indeed, corresponding expan¬ 
sion of the “vertex” operator dr{Q matching the obtained below expansion of the r.h.s. of 
(45) is simple: 

~ + C^di3,k , (46) 

where the index k = 1,..., N in means fc-th branch of Riemann surface (25). 

However, direct substitution of (28) and (46) yields the desirable inhnite set of equations 


{h%o,k=d,G{h\h'^), {h^)p,k=d,G, 


{h\h 


\pV) 


(h* 


]^2,k 


= d,G: 


,{h\h 




only for distinct indices i and fc, because the function G{x,y) has a singularity (see (42)) 
for X = y: 

G{x,y) = Q{x,y) + \n{x-y), (47) 

with the non-singular part 

Qix,y) = ~ + ^a:ln^^ - In (Vl-x^ + ^yl-yA . 

y'l — y^ ^ y — ^ ^ '' 

To complete the construction for i = k, one should observe that substitution of (28) into 
(47) leads to the asymptotic expansion (^ —)■ 0) 

G(ft‘, p<‘)(c)) = Q(h\ +cyf+■■■)+in(pf >+cpf+cyf+...)+in c. 


Since the leading term In ( disappears in (45) after differentiation, we obtain the necessary 
matching expansions for G and dr{(;) for i = k. For instance. 


{h%,k = [Q(h^ h^) + \np[\, {h%ck = I 


(fc) dQ{h^,p) 


dp 


(k)' 
(fc) 

p=h^ 


(k) 

Conservation law densities pm were described in Section 4. Thus, N hrst commuting 
flows are 


(h*)iO,fe 


v/1 - (h^)2 2 h^-1 ^1 - + Vl - {h’^y 


, i ^ k] 


ih%.k 



+ 1 
— 1 


Vi - {h'^y n 

m^k 

n (h^^l - (h”")2 - h'^y/1 - {h’^Y) 

m^k 


X 
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7 Generalized Hodograph Method 


In this Section we finish our procedure for integration of hydrodynamic type system (1), 
also written in equivalent forms (20), (26), (32), (35). 

According to the Generalized Hodograph Method (see detail in [14]), any generic so¬ 
lution r^{x,t) of a semi-Hamiltonian diagonal hydrodynamic type system (5) in a neigh¬ 
borhood of a generic point is given in an implicit form by the algebraic system (4) for 
the unknowns r^{x,t), where ^‘(r) are the velocities of a generic commuting flow (6). In 
arbitrary hydrodynamic variables M*(r) one can easily rewrite the algebraic system (4) (cf. 
[14]) as 

x5l-tVj{u)=w]{u), (48) 

where the hydrodynamic type system (5) has the form 

3 

while commuting hydrodynamic type systems (6) have the form 


3 


In this Section we will modify (4), (48) further in order to get the simplest form suitable 
for the case studied. 

In order to construct solutions of (26) we first need to prove the following result 
suitable for our particular case of system (1) (cf. [11, 12]): 

Lemma: Hydrodynamic type system (26) together with commuting flows (45) has the 
eommon conservation law 


dz 


1 2C'i 

+ J^dt + 


( 1 1 p(h) +1 , pjy) 


dT{r]). 


Proof: If /i —)■ 0, 





p{m -1 


+ 


p(h) 

2(C'0)2(l - 



+/i^(...) + ... 


Taking into account (41) and (44), one can obtain 


1 \ ^ . p{y) + 1 I p{y) 

(^°)V.(,) V2(G0)2 p{r^) - 1 {C^ni-p\p)) 


This conservation law together with (40) can be written in the above potential form, 
where z is a. potential function such that z^ = The Lemma is proved. 

Next we modify the initial Generalized Hodograph formula (4). It should be reduced 
(following the idea in [12]) to a form suitable for our non-diagonal systems (4) and (26). 

Algebraic system (4), (or (48) in arbitrary variables) can be written in the form (here 
di = d/drfl 


X + t 


mo 


diGoo 


(49) 
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where we denoted Ho = (C°)-^Gl = 2C\C^)-^ and 


r 1 I 

2(C'0)2 p(A) - 1 (^0)2(1 - p2(A)) • 

Here we returned to the original parameter A to indicate that now we are working not 
with any particular asymptotic expansion (cf. (28), (44)), but with the original algebraic 
surface as a whole. Indeed, hydrodynamic type system (5) has a conservation law dtHo = 
d^Gi, while the commuting hydrodynamic system has the conservation law drHo = dxGoo- 
This means that ^ diHo ■ rl = 9iGi ■ r* and ^ diHo ■ ^iGoo ■ Taking into 

account (5), (6), (4) and splitting w.r.t. r*, one obtains (49). Multiplying (49) by diHo-dr^ 
and summing up, one arrives at 

xdHo{r) + tdGi(r) = dGooi^)- 

Now we rewrite this equation after the invertible point transformation (r) —)■ (b) as 
xdHo(h) + tdGi{h) = dGooiJo), so the algebraic system (4) becomes 

OHq dGi _ dGao 

^ (96* 56* 56* 

Taking into account diHo = —2(C°)“^, 5jGi = 2(C'°)“^6* — 4C'^(C°)“^ (see (33), here 
5j = 5/56*), we obtain the algebraic system 

^+t(2c'-c%-) = In (yTTy-,)-,yTTy+c»i^ - if?) ■ 

(50) 

where g(b. A) is the inverse function to the function A(b,g) (25). These N equations 
are nothing but the diagonal part of the matrix algebraic system (48). All off-diagonal 
equations are compatible with the diagonal part ([14]). 

So we proved: 


Theorem 1 Hydrodynamic type system (26) has infinitely many particular solutions b^{x, t) 
in the implicit form given by (50) with a free parameter A. 


The algebraic system (50) determines one parametric family of solutions 6*(x, f. A) 
in implicit form and simultaneously g{x,t,X) = —C°(b(x, t. A)). Thus we found one 
parametric family of Hamilton’s equations (8), which are Liouville integrable. 

In fact, using the Generalized Hodograph Method and the nonlinear superposition 
principle implied by this method (see below) we easily obtain multiparametric families of 
integrable metrics. Namely expanding the generating function p{h, A) at different points 
on the Riemann surface p = p{h, A) with the parameters {p, A) (for example when p —>■ — 1 
or p —)■ 6*), one can construct inhnite multiparametric series of new solutions g{h{x,t)). 
Now we demonstrate this idea. To avoid large repeated expressions, we introduce functions 


hR(b, A) = In - g) - + G° 

^ \m=l^ 


Nq \ 

1 + j 
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1. Kruskal series. Substitution of asymptotic expansion (43) into (50) leads to 


Wiih.ii) = In/i + -Pb'/"^^(b) + Pb-f ^(b) + + ..., 

/X 

where, for instance, 

W^~^\h) = C%^ - 2C^. 

ivf>(b) = {Vf - - logC» - log2. 

2. N principal series. Substitution of asymptotic series (27) into (50) leads to 
vy'=>(b, A) = M'.‘,'=>(b) + AH',?’(b) + 
where, for instance, 

= In (^1 + (6^)2 - + {C^6f - b^)^/lTWy. 

Once we found all these coefficients and the Kruskal series hFj^(b), we can con¬ 

struct inhnitely many particular solutions parameterized by arbitrary number of constants 

N oo 

x + t{2C'-C%‘) = J2'E^TWL{h), ( 51 ) 

k=lm=0 

or a functional parameter <yc(A): 

a; + t(2C^ - C%^) = j ip{\)W,{h, X)d\, (52) 

where </?(A) and the contour can be chosen in many special forms. Formulae (51) and 
(52) present the nonlinear superposition principle implied by the Generalized Hodograph 
Method. 


8 Conclusion 


In this paper we considered integrability of semi-Hamiltonian hydrodynamic type system 
(1). We found and presented N inhnite (principal) series of conservation laws and com¬ 
muting flows. This means that one can extract inhnitely many corresponding solutions 
by the Generalized Hodograph Method. Here we did not investigated a completeness of 
conservation law densities (and correspondingly commuting hows). This problem should 
be investigated elsewhere. Here we just mention that (see (30)) for iV = 2n -|- 1, n > 0 


2n+l 

Epf’ 


k=l 


(1 + (6^)2)«-i 
n {b^ - b^) 

m^k 
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This means that corresponding conservation law densities are not linearly independent 
for odd N. Thus one probably should construct additional conservation law densities to 
get the complete set. 

Semi-Hamiltonian hydrodynamic type system (1) holds for the coefficients a^{x,t) of 
a first (polynomial) integral (2). In this paper we found more appropriate set of held 
variables 6^. This natural choice of unknown functions follows also from factorization of 
polynomial expression for the hrst integral with respect to ratio of both momenta®: 


^= 


N N 


n (9 - ■ 


k=0 ^ ' k=l 

Also this hrst integral can be written in two other equivalent forms (see (33) and (36)): 


F=(p2)"exp 


r/k-\-l 


k=0 


{P2f + 


ElV 

B^J 


N/2 


A(g,B), 


where 


q = 


B% 

Pi 


These two equivalent representations for the hrst integral are of big interest: if one hnds 
another hnite-dimensional parameterization of the moments i?^(b) compatible with (34), 
the hrst integral associated with this solution will be no longer polynomial. This more 
general problem should be considered in a separate publication. 


Acknowledgements 

M.V. Pavlov’s work was partially supported by the grant of Presidium of RAS “Funda¬ 
mental Problems of Nonlinear Dynamics” and by the RFBR grant 15-01-01671-a. S.P. 
Tsarev acknowledges partial hnancial support from the grant from Russian Ministry of 
Education and Science to Siberian Federal University (contract N° 1.1462.2014/K). 


References 

[1] M. Bialy, A. Mironov, Rich quasi-linear system for integrable geodesic hows on 2- 
torus, Discrete and Continuous Dynamical Systems - Series A. 2011. V. 29. N. 1. P. 
81-90. M. Bialy, A. Mironov, Cubic and quartic integrals for geodesic how on 2-torus 
via system of hydrodynamic type, Nonlinearity. 2011. V. 24. P. 3541-3554. M. Bialy, 
A. Mironov, Integrable geodesic hows on 2-torus: Formal solutions and variational 
principle. Journal of Geometry and Physics, Vol. 87, No. 1 (2015), P. 39-47. 

[2] Bolsinov, A. V., Matveev, V. S., Fomenko, A. T. Two-dimensional Riemannian met¬ 
rics with integrable geodesic hows. Local and global geometry. Sbornik: Mathematics, 
189(10), 1441-1466, 1998. 

®we remind that = 1 and = g. 


20 








[3] Bolsinov A. V., Fomenko A. T. Integrable Hamiltonian systems: geometry, topology, 
classification. - CRC Press, 2004. 

[4] B.A. Dubrovin, S.P. Novikov, Hamiltonian formalism of one-dimensional systems of 
hydrodynamic type and the Bogolyubov-Whitham averaging method, Soviet Math. 
Dokl., 27 (1983) 665-669. B.A. Dubrovin, S.P. Novikov, Hydrodynamics of weakly 
deformed soliton lattices. Differential geometry and Hamiltonian theory, Russian 
Math. Surveys, 44 No. 6 (1989) 35-124. 

[5] M. V. Deryabin, Polynomial integrals of dynamical systems and the Lax reduction. 
Math. Notes, 61 No. 3 (1997), 363-365. 

[6] J. Gibbons, S.P. Tsarev, Reductions of the Benney equations, Phys. Lett. A 211 
No. 1 (1996) 19-24. J. Gibbons, S.P. Tsarev, Conformal maps and reductions of the 
Benney equations, Phys. Lett. A 258 No. 4-6 (1999) 263-271. 

[7] V. V. Kozlov, Polynomial integrals of dynamical systems with one-and-a-half degrees 
of freedom. (Russian) Mat. Zametki 45 No. 4 (1989) 46-52; translation in Math. 
Notes 45 No. 4 (1989) 296-300. 

[8] M. V. Pavlov, Algebro-geometric approach in the theory of integrable hydrodynamic 
type systems. Comm. Math. Phys. 272 No. 2 (2007) 469-505. 

[9] M. V. Pavlov, Classihcation of integrable hydrodynamic chains and generating func¬ 
tions of conservation laws, J. Phys. A: Math. Gen., 39 No. 34 (2006) 10803-10819. 

[10] M. V. Pavlov, Integrable hydrodynamic chains associated with Dorfman Poisson 
brackets. arXiv: 1008.4530. 

[11] M.V. Pavlov, S.P. Tsarev, Tri-Hamiltonian structures of Egorov systems of hydro- 
dynamic type, Funct. Anal. Appl., 37 No. 1 (2003) 32-45. 

[12] M.V. Pavlov, S.P. Tsarev, Classical Mechanical Systems with one-and-a-half De¬ 
grees of Freedom and Vlasov Kinetic Equation, in: Topology, Geometry, Integrable 
Systems, and Mathematical Physics, AMS Translations: Series 2, vol. 234, 2014, pp. 
337-371. 

[13] Rozhdestvensky B. L., Yanenko N. N. Systems of quasilinear equations and their 
applications to gas dynamics. Translations of Mathematical Monographs, 55. AMS, 
Providence, RI, 1983. 

[14] S.P. Tsarev, On Poisson brackets and one-dimensional Hamiltonian systems of hy¬ 
drodynamic type, Soviet Math. Dokl., 31 (1985) 488-491. S.P. Tsarev, The geometry 
of Hamiltonian systems of hydrodynamic type. The generalized hodograph method. 
Math. USSR Izvestiya, 37 No. 2 (1991) 397-419. 


21 


